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.'^ . Abstract 

^^ ' According to extensive ab initio calculations of lattice QCD, the very large energy 

^^ , density available in heavy-ion collisions at SPS and now at RHIC must be sufficient 

to generate quark-gluon plasma (QGP), a new state of matter in the form of plasma 

O of free quarks and gluons. The new state of matter discovered at RHIC seems to be 

.^ ' perfect fluid rather than free plasma. Its shear viscosity is assumed to be almost zero. 

In this work, I first considered the theoretical and phenomenological consequences of 

this discovery and finally asked questions about the nature of phase transition and 

properties of matter. It is important to answer these questions, otherwise QGP will 

^ ' remain a kind of black box; one sends a signal via new experiments or simulations 

[^ , or models and gets another one from it. I will show that some promising ideas have 

Cn ■ already been suggested a long time ago. I will also suggest a new phase diagram 

^-P , with separated deconfinement and freeze-out boundaries and a mixed state of thermal 

quark matter and bubbles of hadron gas. 
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1 Introduction 



Theory that describes the dynamics of phase transition to quark-gluon plasma (QGP) 
and back to hadronic matter is still failed. In high-energy experiments, we can only 

,^ . measure produced particles in the final state, i.e., after chemical and thermal freeze-outs. 

Consequently, occurrence of equilibration processes while energy density decreases can 
not be confirmed experimentally. Same statement is valid for degrees of freedom during 
hadronization processes. It is known that thermal models work well in the final state. They 

5^ \ presume a charge-conserved hadronic phase and have no access to phase transition. In 

other words, the phenomenology of QGP and particularly the dynamics of phase transition 
seems to represent a kind of black box. This might explain why we used to assume that the 
deconfinement temperature Tc at small chemical potential is coincident with the chemical 
freeze-out temperature Tfo pp. What are the consequences of this assumption? For 
example, what type of strongly interacting matter we have in mind that undergoes a 
phase transition back to confined hadrons, very much rapidly expands and finally freezes 
out through chemical processes without any change in its temperature? Is it possible to 
describe it by QCD? Or should we assume that the chemical freeze-out does not exist? 

Although the energy density available in heavy-ion collisions at CERN-SPS T and now 
at BNL-RHIC exceeds the critical value calculated in lattice QCD (ec ~ 0.7 GeV/fm~^), 
many physicists still debate on QGP signatures under these laboratory conditions. If 
Tc = Tfo should hold at small chemical potentials, phenomenological signals that char- 
acterize phase transition should remain measurable in the final state. At least QGP 
signatures that are not sensitive to a medium, such as color screening or J/^ dissocia- 
tion into two leptons have to survive at Tc = Tjo- On the basis of current experimental 
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and theoretical progress, we can think of other solutions, like Tc >> Tf^- If Tc >> Tj^, 
then QGP signatures are very much contaminated (in-medium-modified) so that the final 
state does not reflect QGP production. Furthermore, we can assume that the production 
of QGP requires much higher temperature than the currently available temperature at 

RHic m. 

QCD predicts that matter under extreme conditions is simple but it does not say 
anything about matter being ideally formed. For temperatures larger than Aqcd, the 
matter is in a deconfined phase. Only a good understanding of the dynamics of phase 
transition will help us to characterize QGP. 

2 Transition to new phase(s) 

The end of the hadron era was predicted fifty five years ago [5|. No doubt that hadrons 
will go into a new state of matter at sufficiently high temperatures and densities p^ . The 
wide acceptable framework to study the phase transition of strongly interacting matter is 
given by QCD. QCD predicts that the hot hadronic matter becomes simple. This is not 
necessarily weakly interacting. QCD has been studied on a lattice for the past thirty years; 
QCD Lagrangian has to be discretized and everything has to be put on a finite a space- 
time lattice. It has been shown that a change in the phase of matter undoubtedly exists at 
sufficiently high energy densities. The degree of freedom markedly increases in a relative 
narrow region of temperatures. For a purely gluonic system, for which the equation of 
state can be computed without approximations, there is a deconfinement phase transition. 
It is of the first order and the critical temperature is T^ ~ 270 MeV [7j. 

The difficulties in lattice calculations start when the fermion sector is switched on. 
Nevertheless, it has also been observed that the chiral symmetry is restored at the same 
critical temperature Tc ~ 154 — 174 MeV (depending of quark flavors) as that of decon- 
finement transition. The restoration of chiral symmetry means that the effective mass of 
quarks forming the hadron states becomes zero. Another important consequence of chiral 
symmetry breaking restoration is the disappearance of the mass degeneracy of hadronic 
states with the same spin but different parity quantum numbers. Dynamical quarks can 
only be included in lattice QCD in certain approximations. The order of phase transition 
in full QCD is not yet completely known. We merely see a rapid change in bulk thermo- 
dynamic quantities. Consequently, the transition is known as cross-over. The behavior of 
QCD phase transition at a finite chemical potential is not yet known from first principle 
of QCD. Nevertheless, if we look at the lattice results at a zero chemical potential. Fig. ^ 
we find that the system at temperatures of 4 — 5Tc remains below the Boltzmann limit [Hj; 
^SB ~ gir'^T^/SO. It seems that e/T^ will remain constant at higher temperatures. This 
means that the deconfined system is still strongly correlated. This is also valid at finite 
chemical potentials. 

Apart from the known properties of lattice calculations, like heavy quark masses and 
finite lattice size and spacing, we have from this figure evidence that hadrons undergo 
a transition to a strongly correlated phase. Another important finding is the nature of 
phase transition. It is not a real phase transition but a rapid change in the degree of 
freedom. Characterizing the nature of matter above Tc is essential for particle physics. It 
is important in order to understand the development of the universe in early stages. 
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Fig. 1: Energy density in lattice QCD calculated for different quark flavors. It is clear 
that the energy density at the temperature 1.2Tc < T > 4Tc is smaller than the energy 
density in the Stefan-Boltzmann limit. 



3 RHIC results and lattice response 

Using elliptic flow, it was possible to conclude that RHIC has produced thermalized matter 
at a very high energy density ^. For the first time, hydrodynamics with a zero viscosity 
can describe heavy-ion reactions. In relativistic hydrodynamics, under the assumption of 
local thermalization, the number density, current and energy-momentum tensors are 
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where r] and C are the shear and bulk viscosities, respectively. The number density and 
entropy density are conserved. In strongly interacting quantum fields, the ratio of shear 
viscosity to entropy density is as |3] 
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As mentioned above, it was possible to describe RHIC results by hydrodynamics, even if 
the last two terms in Eq. ^are entirely removed. This means that matter above Tc should 
also have a zero viscosity. It is a fluid rather than a free plasma; the correlations do not 
completely vanish. 

The pioneering quenched lattice calculations performed by Atsushi Nakamura of Hi- 
roshima University for the transport coefficients of gluonic plasma ^U] seem to confirm 
the above results. The lattice calculations are performed in the framework of the linear 
response model. 
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^The official press release had the title: "RHIC scientists serve up perfect liquid". This mainly has been 
based on a comparison with hydrodynamics. Because in the hydrodynamic codes, one has to introduce 
"many" inputs, e.g., initial conditions, equation of state, phase transition, chemical freeze-out and thermal 
local equilibration in order to solve the equations of motion numerically and describe the hydrodynamical 
evolution, one has to be very careful with this kind of comparison. 



< Ti2{x,t)Ti2{x',t') >ret IS related to Green's function of the energy-momentum tensor. 
Fig. [21 depicts the results for temperatures up to 30Tc. Shear viscosity in the perturbation 
theory is as jTJ 
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where r/i = 27.126, fj,* /T = 2.765, fi = Tc/Aqcd and g is the running coupling constant. 
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Fig. 2: Top: rj/s obtained by quenched lattice simulations (symbols) compared with 
perturbative results (lines). Bottom: shear viscosity in GeV^ units, /u is a parameter 
given in Eq. 21 

As mentioned above, only the gluonic sector has been taken into account in these 
lattice calculations. The fermionic part entirely vanishes. In this case, the bulk viscosity 
is almost zero. It is clear that lattice calculations agree with the perturbation theory at 
high temperatures. However, for T < 3Tc, we find that the lattice results lie below the 
perturbative lines. 

It is needed to include quarks. By doing this, we can describe the transport properties 
of matter above Tc and make a better comparison with phenomenological results. Also, 
we need to improve the inputs for hydrodynamical codes. 

The collision energy is sufficiently high to catapult hadrons into a new phase. Thus, we 
now have the situation that the properties of recently discovered new state of matter has 
to be confirmed by three independent methods: analytical, numerical and experimental. 
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Fig. 3: J /'^ bound state has been found to survive at almost twofold Tc- Curves cited 
from .12.. 



QCD-like models and/or effective models can be used in the analytical method. To 
the author's knowledge, the only analytical method applied so far - apart from limitations 
- is in [^ . 

The most powerful numerical method is lattice QCD with dynamical quarks, which still 
uses nonphysically heavy or almost vanishing quark masses and faces a serious challenge 
when the chemical potential becomes finite. As mentioned above, the properties of new 
state of matter are studied in lattice QCD in quenched approximation jJOl- The lattice 
results have to be interpreted with respect to all these approximations. For example, the 
critical energy density and temperature at the physical quark masses are larger than those 
given above, namely, T^ — 200 MeV, ec — 2 GeV/fm~^. This has the consequence that 
Tc 7^ Tfo. The deconfinement point is not the same as that of chemical freeze-out ll|. 

The experimental method is the ultimate goal. The proof that a new phase of matter 
is produced - beyond any doubt - has to be adduced experimentally. Experimentalists 
have maintained this position since more than two decades. Particle theorists used to 
tell them that hadrons for T > Tc will go into an "ideal" deconfined phase and into free 
plasma. 

Plasma is a phase of matter in which charges are screened due to the existence of 
other mobile charges. This will modify Coulomb's Law. Recent lattice results on spectral 
functions show that the J/^ bound state can survive even at ~ 2Tc, (Fig. EJ. In these 
quenched lattice calculations, Tc — 270 MeV. As mentioned above, RHIC results imply 
that the new state of matter is a fiuid rather than a free gas. 

On the other hand, when we think about the situation that, in heavy-ion experiments, 
one starts and ends up with hadrons, we realize that the experimental proof is non trivial. 
Hadrons will be accelerated, collide against each other and freeze out to produce new 
"hadrons". The latter will be counted by detectors. This should illustrate the importance 



of making precise theoretical predictions. Nevertheless, I think that experimental physics 
is now studied using powerful tools. 



4 Gluon condensates and correlations 
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Fig. 4: Gluon condensate of SU{3) (solid line) and an ideal gluon gas (broken line) 
compared with the cases where light (heavy) dynamical quarks are included, denoted by 
open (full) symbols, respectively. Curves cited from [T? 



Many ideas about the complex structure of phase transition have been suggested. For 
instance, the soft modes of confined hadrons can survive the restoration of chiral symmetry 
breaking [T3j and that hadron plasma can be formed instead of QGP [14] . Also there are 
many lines of evidence that deconfined matter don not behave like free gas. For example, 
gluon condensates |151ll6j remain finite above Tc- Gluon condensates have been calculated 
at a finite temperature and in the presence of massive quarks ^HI- It has been assumed 
that such quarks due to strong interactions will modify the thermal properties of gluon 
condensates. Therefore gluon condensates remain almost unchanged with the half of their 
value in the confined phase ^Sj. This means that matter seems to be sticky J17j . 

(G2)t = {G\ + mg{i;g4^g)o - m,{i,,,Pq)T - {%)t. (5) 

Fig. ^depicts this quantity. At T > Tc, the condensates at a zero chemical potential seem 
to be T— independent. 

Thus, we have concluded that hadronic matter is expected to undergo phase transition 
to a new state of matter at high temperatures. The new state of matter is strongly 
correlated due to the non vanishing gluon condensates and the existence of finite hadronic 
modes. 

5 Phase space and QCD phase diagram 

Phase transition at a very large chemical potential is most likely of the first order; at Tc, 
only one phase can exist. This first-order line is expected to end up with a critical end 



point of the second order. Its location was the subject of different lattice simulations. To 
date, there has been no final result. In dynamical QCD with 2 + 1 flavors of staggered 
quarks of physical masses ^H]) the endpoint has been localized at T = 162 it 2 MeV and 
^b = 360 lb 40 MeV. According to the same reference the critical temperature at ^f, = 
is 164 lb 2 MeV. From recent lattice simulations UH] and comparisons with the resonance 
gas model PHI, we assume that Tc{fib = 0) should be ~ 200 MeV. The temperature of the 
chemical freeze-out at a zero chemical potential is ~ 174 MeV. Another determination of 
the endpoint has been reported in j21j : ^^ ~ 420 MeV. The corresponding temperature 
has not been estimated. 

We can summarize that the endpoint is at approximately ;Uf, ~ 400 MeV or fJ-b/Tc ~ 2. 
This value corresponds to collision energy y/SNN ~ 9 GeV. The lead beam at 40 AGeV 
(SPS) can be used to scan this region. Note that, at this energy, we have fully unexplained 
peaks in the ratios of strangeness to non-strangeness particle yields |22| . That the ratios 
of strangeness hyperons to non-strangeness hadrons (pions) all have a maximum value at 
the same energy is an indication of strangeness asymmetry. This can only be achieved in 
the plasma phase |5^. i.e., it assumes that the interacting system at this energy should 
undergo phase transition to QGP. The phase transition is likely of the first order. 

On the other hand, we do not observe such a peak at higher energies. It is an open 
question, why QGP-signatures e.g., color screening and strangeness asymmetry are not 
seen, when colliding energy increases? It is believed that, at RHIC and SPS |21] energies, 
we have produced the highest thermalized system ever known. The reason might be 
the nature of system produced and the dynamics controlling phase transition. In other 
words, the new state of matter might not have the properties of QGP, for which we have 
suggested phenomenological signatures. 



One possibility to explain the difficulties with detecting phenomenological signatures 
at low chemical potential is the non-equilibrium quark occupancy of phase space, i.e., 
7j 7^ 1. At finite temperature T, strangeness fj-s and iso-spin /i/3 and baryo-chemical 
potential fiB, the pressure of one hadron is 
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where e{k) = {k'^ + m?)^''^ is single-particle energy and ib stands for bosons and fermions, 
respectively, g is spin-isospin degeneracy factor. A = exp(/x/T) is fugacity parameters. 
jjL is the chemical potential multiplied by corresponding change. The total pressure is 
obtained by summation all hadron resonances. 7 appears in front of Boltzmann exponen- 
tial, exp(— e/T). It gives averaged occupancy of phase space relative to equilibrium limit. 
Therefore, in equilibrium limit, 7 = 1. Assuming time evolution of system, we can describe 
7i as ratio between the change in particle number before and after chemical freeze-out, i.e. 
li = ni{t) / ni{co) . The chemical freeze-out is defined as time scale, at which there is no 
longer particle production and the collisions is entirely elastic. In case of phase transition, 
7j is expected to be larger than one, because of non-equilibrium processes, large degrees 
of freedom and expanding phase space. 

The dependence of single-particle entropy on collision energy is related to averaged 
phase space density. In Boltzmann limit and for one particle, we get that 
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where e is single-particle energy. In this expression, s/n directly relates to e. s is the 
entropy density, n is the particle number density. Apparently, s/n gets a maximum value. 
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Fig. 5: Entropy per particle s/n as function of ^snn- Singularity is located at almost 
the same energy as the peaks of particle ratios of strangeness to non-strangeness |22) . 

when ^ = e. Because of Boltzmann limit, the maximum, in this classical system, is unity. 
Depending on n, we can insert particles into phase space. Maximum occupation is reached 
at /x = e. Beyond this limit, it is prohibited to insert more particles. On the other hand, we 
can expect - at least theoretically - that the occupation value is larger than one, if phase 
space itself is changed. The latter situation is most likely provoked by phase transition. 

The connection between this theoretical discussion and particle ratios of strangeness 
hyperons to non-strangeness is given by 7j. As in 22 , s/n plays the same role as 7j. 

The results of s/n vs. ^/snn are depicted in Fig. El Full grand-canonical statistical set 
of the thermodynamic parameters is used. In this case, the complete dependence of s/n 
on T and fi and consequently on ^s^n, can straightforwardly be obtained by deriving s 
and n from Eq. El 

For 7g = 7s = 1, we find that s/n increases as the energy raises from AGS until 
low SPS energies {^snn < 9 GeV). For higher energies, s/n remains constant. This 
behavior might be an indication to strong compensation of collision energy in this region. 
Although, more energy is introduced to system, number of particles allowed to occupy 
phase space remains constant. It is an indirect signature that the phase space itself 
remains constant. 

For varying 7^ and 7^, we find a singularity at ^snn ~ 7 GeV. Equivalently, phase 
space is maximum at this energy. At higher energies, s/n decreases. Although, energy 
increases and consequently the produced particles, single-particle entropy decreases. 
This means that phase space shrinks. At RHIC energies, the shrinking becomes slower 
than at SPS. If this model gives the correct description, we now might have for the 
first time a theoretical explanation for dependence of phase space on energy. The phase 
space at SPS energy is apparently larger than at RHIC. Same behavior has been found 
experimentally ^17^. The consequences are that QGP might be produced at SPS and 
detecting its signatures at RHIC might be non trivial. 



The situation is different at small chemical potentials. It is most likely that the tran- 
sition can not be characterized by a line. As shown above, the transition in this region is 
a rapid cross-over. The critical temperature is an average, (Fig. P). Therefore, we expect 



Quark-Gluoii Plasma (QGP) 

phase ^ 

^G+Qg\ < t^j/j >- 

^endpoint of second order 




Quark- Matter 

first-order phase raiisition 



l^b[GeV] 



Fig. 6: Phase diagram of QCD. The values of the chiral symmetry order parameter < 
■ipip > are also given. We distinguish between the freeze-out and transition lines. At small 
chemical potentials, the hadrons most likely undergo a smooth transition to a thermal 
phase accommodating a hadron gas and a fluid of quarks and gluons. At liight T, the 
mixed phase entirely disappears. Matter will then be formed in QGP. 



that the transition is a kind of wide boundary. From Fig. ^ we can roughly estimate that 
the transition boundary has a width of at least ATc- At temperatures up to 4Tc, the degrees 
of freedom in the new phase of matter are 30% smaller than that in the Stefan-Boltzmann 
limit, (Fig.d). 

There are indications that the hadronic bound states can survive above Tc, (Fig. ^. 
We have discussed the hydrodynamical results at RHIC energies. From all these results, 
it seems physically consistent to suggest that the new phase of matter is a mixture of 
the surviving hadron gas and a fluid of quarks and gluons [25j . The hadron gas might 
form bubbles inside the quark fluid. This explains two important features, the degrees of 
freedom and existing hadronic bound states above Tc. 

Based on the above discussion, we suggest a "new" phase diagram in Fig. El We 
separate the line of equilibrium chemical freeze out from that of phase transition. At 
a zero chemical potential, phase transition takes place at T ~ 200 MeV ^ E]- The 
equilibrium chemical freeze-out takes place at T ~ 174 MeV |lj. The two lines slightly 
become shorten with increasing chemical potential. Phase transition at chemical potentials 
up to the value corresponding to the endpoint is no longer a cross-over. Consequently, 
QGP is expected at much higher temperatures than Tc- Between the hadron gas and 
the QGP there is a mixed phase, in which the two degrees of freedom partly exist. The 
properties of matter in the mixed phase are as follows. 

It is an ordinary "nearly perfect" fluid (it is not free gas or plasma). 

It has very small but finite viscosity (it is not a superfluid). 

It has electric resistance (it is not a superconductor). 

It is a strongly correlated matter, since hadronic bound states can survive above T^. 

It is a sticky matter, since gluon condensates remain finite above Tc- 



In the light of these properties, we need to answer the following questions: 

• What is the order of phase transition at small chemical potentials (related to RHIC, 
LHC or early stage of universe)? What is the order of phase transition from the 
mixed phase to QGP at very high temperatures? 

• How can we confirm the transition at small chemical potentials phenomenologically? 

• What are the in- medium modifications of hadron properties at finite T and Hb j2f)j ? 

• What are transport properties above Tc? 

We have to answer all these questions before we face a much difficult challenge with CERN- 
LHC and GSI-FAIR. We have to have clear theoretical descriptions and non ambiguous 
phenomenological signatures describing the QCD phase diagram. Non-perturbative QCD 
is a good candidate to answer some of these questions. Phenomenological studies are also 
very essential. 
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